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POST BUCKLING ANALYSIS OF SANDWICH BEAMS WITH FUNCTIONALLY GRADED 
FACES USING A CONSISTENT HIGHER ORDER THEORY 

AKIL AHMED 



In this project work, a new consistent higher order theory has been developed for the buckling and post buckling 
analysis of sandwich beams having functionally graded faces. The model combines an approximation for the displacement 
fields with third order in the thickness direction. The governing equilibrium equations of motion and the related boundary 
conditions have been derived. The parametric study of the beam behaviors under different situations is performed and the 
important influencing parameters are identified. 
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Functionally graded materials (FGMs), a new generation of composites proposed in the earlier 1980s are achieved 
by two or more constituent phases with a continuously variable composition along the thickness direction. It has potential 
application in different field of engineering, especially in defense and aerospace. FGMs have drawn considerable attention 
as they are made from a mixture of metals and ceramics and are further characterized by a smooth and continuous change of 
the mechanical properties from one surface to another [1]. Some research has been carried out to analyze the buckling and 
post-buckling response of structures made of FGMs. But more effort in the research is required for more accurate analysis to 
explore properties and response under different types of environment and different type of configuration. Yaghoobi and 
Torabi [1] studied post-buckling and nonlinear vibration analysis of geometrically imperfect beams made of FGMs resting 
on nonlinear elastic foundation subjected to axial force. 

Loja et al [2] studied the static and free vibration behavior of functionally graded sandwich plate type structures, 
using B-spline finite strip element models based on different shear deformation theories. The effective properties of 
functionally graded materials are estimated according to Mori-Tanaka homogenization scheme. Simsek and Reddy [3] 
proposed a unified higher order beam theory which contains various beam theories as special cases for buckling of a 
functionally graded (FG) micro beam embedded in elastic Pasternak medium based on the modified couple stress theory 
(MCST). This non-classical micro beam model incorporates the material length scale parameter which can capture the size 
effect. The non-classical beam model reduces to the classical beam model when the material length scale parameter is set to 
zero. Li and Batra [4] derived analytical relations between the critical buckling load of a FGM Timoshenko beam and that of 
the corresponding homogeneous Euler-Bernoulli beam subjected to axial compressive load for clamped-clamped (C-C), 
simply supported-simply supported (S-S) and clamped-free (C-F) edges. Simsek et al [5] developed a micro scale 
functionally graded Timoshenko beam model for the static bending analysis based on MCST. They assumed the material 
properties of the FG micro beams to vary in the thickness direction and are estimated through the Mori-Tanaka 
homogenization technique and the classical rule of mixture. 
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Simsek [6] performed non-linear dynamic analysis of a FG beam with pinned-pinned supports due to a moving 
harmonic load by using Timoshenko beam theory with the von-Karman's non-linear strain-displacement relationships. 
Huang and Li [7] performed an analysis of the stability of circular cylindrical columns/beams composed of FGM is made 
where shear deformation is taken into account. The coupled governing equations for the deflection and rotation are derived 
based on the traction-free surface condition and a single governing equation is further obtained. Singh and Li [8] presented a 
mathematical model to compute the buckling loads of uniform and non-uniform functionally graded columns in the axial 
direction. The columns with spatial variation of flexural stiffness, due to material grading and/or non-uniform shape, are 
approximated by an equivalent column with piecewise constant geometrical and material properties. Kang and Li [9] 
investigated the mechanical behaviors of a non-linear FGM cantilever beam subjected to an end force are investigated by 
using large and small deformation theories. Young's modulus is assumed to be depth-dependent. Yang and Chen [10] 
presented a theoretical investigation in free vibration and elastic buckling of beams made of FGMs containing open edge 
cracks by using Bernoulli-Euler beam theory and the rotational spring model. They assumed that the material properties 
vary along the beam thickness only according to exponential distributions. 

METHODOLOGY 

A new efficient model has been developed for the buckling and post buckling analysis of sandwich beams 
consisting of functionally graded faces. The model approximates for the axial displacement fields with third order 
assumption. The model approximates for the materials properties in the thickness direction. The governing equilibrium 
equations of motion have been derived by using Hamilton's Principle. Analytical solutions are obtained for the beam for 
buckling and post buckling response under the loads. Parametric study is performed for the present theory 

Consider a sandwich beam with functionally graded faces as shown in the figure (Figure 1) having length L, 
width b, and total constant thickness h. The following distances are taken from the mid surface of the beam: 

Distance of bottom fibre of the bottom face = hi 

Distance of bottom fibre of the core = h 2 

Distance of top fibre of the core = h 3 

Distance of top fibre of the top face = h 4 
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Figure 1: Sandwich Beam with FGM Faces 





The constituent materials are taken to be ceramic and metal and its volume fractions of the ceramic and metal 
follows the power law [11] as: 

V c = ( z hl \ for bottom face 

\h 2 — hi / 
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V c — 1 for core 

_ ( z ^ \ r Qr fo ott:om f ace 

c \h 3 -h 2 ) ' ' 

V m = l-V c (1) 

The subscripts m and c in the above equation mean to the metal and ceramic constituents, respectively and z is the 
coordinate in the thickness direction of the beam. The superscript k is the volume fraction index that varies greater than or 
equal to zero. For k = 1, the variation is linear for both ceramic and metal and for k = 0, the beam is fully ceramic. 

The mechanical properties of FGM beam here are determined as per the volume fraction of the ceramic and metal 
assumed from the power law. The modulus of elasticity E and density of material p assumed in the thickness direction z 
based on Voigt's rule [11]; and Poisson's ratio v is assumed to be constant [12] as expressed here 

E = E{z) = E C V C + E n (l- V c ) 

P = P(z) = P C V C + P m (1 - V c ) 

v (z ) = v„ (2) 
The displacement field is assumed from the Reddy's third order shear deformation theory [13] 
u(x, z) =u 0 +z [0 O - fcz 2 (0 o + 
w(x, z) — w 0 

k = A (3) 

3h 2 v ' 

In the above equations, u and w are the displacement components in the x and z directions; cp 0 is the rotation of the 
transverse normal about y axis; the generalized displacements are u 0 , w 0 , <Po- 

The strain displacement relations are assumed as follows 

du 1 ,dw. 2 

£x ~J^ + 2^ 

dw . du ... 

Vv 7 = — + — (4) 

xz dx dz v ' 

The constitutive equations are taken in the formulation as following 
Ox = Qn(z)£x 

Tzx = Qss(z)Yxz (5) 
Where 

E(z) 



(l-v 2 ) 



&sW = 5^ (6) 

The potential energy (P) and work done (W) [14] by the load can be written as 
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U = \ S v iPx e *z + t«v„) dv (7) 
W = \i L 0 P^fdx (8) 
Now, the equilibrium equations are derived with the help the Hamilton's Principle [15] 

8 J*(T - U + W) dt = 0 (9) 
The equilibrium equations are derived as follows: 
2Pw xx + 6bk[D 55 0 >x - 3kF 55 ( 0 rX + w xx )] + 
2b[A u u xxx + B ±1 0 xxx - kE^{0 xxx + w xxxx )] = 0 
An u xx + Bn0,^ - fc£»n( 0, xx + w jXXX ) = 0 

,4 55 0 - 3fc£> 55 (0 + w x ) + 9fe 2 F 55 (0+ w, x )=0 (10) 
Where 

f [^+(/l4-/l3] 2 3 

-[— +(hi-h2] 

.^+(/l4-/l3] - . 

^55.055.^55 = /V.,,, h7l (l,Z 2 ,Z 4 )Q 55 dz (11) 
-[— +(ftl-fl2] 

The following displacement fields [16] are assumed for the first buckling mode for simply supported beam: 
w(x) = a sin (^) 

0(x) = 0 cos (f) (12) 
For clamped-clamped boundary condition, the following displacement [18] field is taken for analysis: 
w(x) = cos f^p) - 1 

0(x) = scos(^)- 1 (13) 
For simply supported beam, the critical buckling load calculated as 

n ^ , ,n ^ 2 ,n,,2^ 1 r (9fc 2 F55-3K-D 55 ){(l86fc 2 F 55 -66feD 55 )+^(26fcD 11 -26fcE 11 )L 

P = (bkDii - bkEn) — + 9bk A F 55 - [ = — ] 

L V 11 117 i. 2 55 J L /l 55 -3fcD 55 + 9fc 2 F 55 J 

RESULTS 



We analyze a sandwich beam with functionally graded faces having alumina (ceramic) and aluminium (metal) as 
shown in the figure (Figure 1) with Ec = 380 GPa, Em = 70 GPa, v = 0.23 [17]. The top and bottom faces of the sandwich 
beam are taken 0. 1 h each. The analysis is performed on critical buckling load with the effect of aspect ratio of the beam on 
it and post buckling response. It is observed that shear plays a major role on critical buckling load and post buckling 
response as well. In Figure 2, the post L/h versus critical buckling load is plotted for simply supported beam; 
whereas variation of buckling load is plotted against length to thickness ratio of for clamped-clamped beam in Figure 3. 
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Figure 2: Variation of Buckling Load with L/h for Simply Supported Beam 




Figure 3: Variation of Buckling Load with L/h for Clamped-Clamed Beam 

CONCLUSIONS 

A consistent higher order theory for the analysis of sandwich beams having functionally graded faces has been 
developed here. The theory is formulated and assessed for different parameters. These parameters have significant 
influences on the mechanical behaviour of this type of structures. Since the face sheets made of FGM are suitable for the 
structures often subjected to the rapid variation of temperature, the buckling and post-buckling analysis is very important 
for analysis and design of this type of structures. Analytical solutions for the beam with functionally graded faces will be 
useful to assess the theory. This study will be of an important step in the development of functionally graded structures. 
The post-buckling behavior of a sandwich beam with functionally graded faces is investigated based on a higher-order 
shear deformation theory. The governing equations and the boundary conditions are derived using the Hamilton's principle 
and the governing equations are solved by closed-form method. The effects of length-to-thickness ratios and boundary 
conditions on buckling and post-buckling behaviors are analyzed in this study. From the analysis, it is observed that shear 
has the major effect on the critical buckling load and post-buckling response For the thin beam having high length-to- 
thickness ratios, the effect of shear deformation is negligible. In addition, the clamped-clamped beams take more axial load 
compare to the simply supported beams. 
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